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Abstract 

A relaxed factorization is used to obtain many of the properties obeyed by the confluent 
hypergeometric functions. Their implications on the analytical solutions of some interesting 
physical problems are also studied. It is quite remarkable that, although these properties 
appear frequently in solving the Schrodinger equation, it has been not clear the role they play 
in describing the physical systems. The main objective of this communication is precisely to 
throw some light on the subject. 

1 Introduction 

There is no doubt on the main role played by special functions in theoretical and mathemat- 
ical physics. In general, they are used to simplify the original problem by transforming its 
mathematical description from a rather cumbersome form into a simpler and well known one. 
The physical scenery is thereby clarified: the solutions of the simplified problem can be easily 
analized and their most relevant qualitative features could be depicted in terms of the involved 
parameters. As a matter of fact, a considerable amount of basic research has been developed in 
the study of the differential equations obeyed by special functions, mainly in their connections 
with the problems appearing in all the branches of the physical theories. Of particular interest 
are the hypergeometric and confluent hypergeometric functions (h.f. and c.h.f. respectively), 
in terms of which almost all the solutions of exactly solvable problems in quantum mechan- 
ics can be written (e.g. those related with the linear and p-dimensional harmonic oscillator, 
hydrogen-like, Poschl- Teller, Wood-Saxon, Hulthen, Morse, Eckart and Scarf potentials among 
others) . The standard mechanism in these cases takes into account appropriate transformations 
of the involved variables and functions from the Schrodinger into an hypergeometric or confluent 
hypergeometric equation (h.e. and c.h.e. for the last two respectively) Remark on the fact 
that the mathematical properties of the h.f. and c.h.f. lead then to the characterization of some 
important physical features of the solutions. That is the case, for example, of the quantization of 
the energy eigenvalues by imposing the Schrodinger solutions of classically confined systems as 
square integrable functions. On the other hand, an interesting perspective arises by considering 



the procedure in the reverse order: one departs from the h.e. or the c.h.e. and the conditioning 
transformations leading to the Schrodinger equation are to be determined (see important works 
by Natanzon Q and by Nikiforov and Uvarov ||). Such a procedure allows the identification of 
a wide set of functions that can be understood as physically meaningful potentials. 

Considering all this theoretical richness, it seems apparent the presence of relations con- 
necting the stationary Schrodinger solutions of diverse quantum systems. Although they have 
been frequently mentioned in the literature since long time ago, there is still lacking an updated 
exposition of the topic. Quite recently, the present authors have reported results which could 
shed a new light on this subject [0]. In that work a refined factorization was applied to analize 
the c.h.e. and its solutions, as well as the implications of their mathematical properties on the 
wave- functions of some interesting physical problems. The refined factorization goes deeply into 
the possibilities of the conventional factorization method providing with additional significative 
information ||. 

In this contribution we shall discuss some of the main results published in Q whereas it 
is also reported more on their physical consequences. The Section 2 is devoted to the refined 
factorization of the c.h.e. The action of the involved factorization operators on the confluent 
hypergeometric parameters underlies the analytical and algebraic properties of the c.h.f space. 
Section 3 deals with the mapping from the c.h.e. into diverse Schrodinger equations. The 
identification of the involved conventional physical potentials gives the chance to translate the 
results obtained in Section 2 from a mathematical into a physical language. 



2 Confluent hypergeometric intertwiners 

We shall work on second order differential operators of the form 

L( a ,c) = x-r^ + (c-xj— -a (1) 
v ' ax z ax 

where the pair (a, c) represents a point on R 2 . Once the values of a and c have been given, the 
kernel elements of Lu c \ can be characterized as c.h.f., i.e., f(a, c; x) is a c.h.f. iff Lu c \f(a, c; x) = 
0, or in other words, f(a,c;x) G K-ua- Thereby, one can look for a way to connect the kernel 
K.( a ,c) with another one K-fa,£) an d the conditions relating their respective parameters (a, c) and 
(a, c). With this aim let us consider an arbitrary differential operator X, defined by its action 
on a kernel element f(a, c; x) of Lu c \, as follows 

Xf(a, c; x) := X {a>c) f(a, c; x) = £(o, c; x) (2) 

where, according with our earlier convention, £(a,c;x) is a c.h.f. iff £(a,c;x) S £(a,c) an d the 
pair (a, c) depends on a and c. We shall say that the operator X is in a free index notation 
whereas it is written in the confluent hypergeometric notation Xr a ^ in dependence on the kernel 
K-Uc) on which it is acting. The benefits of this notation become clear in the composition of 
any pair X, Y of these operators. For instance, if Y is such that Yf(a,c;x) = g(a,c;x), then 
the action of XY on the kernel K>( a ,c) reads 

(XY)f(a,c;x) := X (5)g) (Y M f(a,c;x)J = X {h ^g(a, c; x). 
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On the other hand, the action of YX on the same kernel gives 

(YX)f(a,c;x) :=%,<;) ( x (a,c)f( a , c ; x )) = Y (a,c)C(a,c;x). 

Taking full adavantage of our considerations we look now for a couple of differential operators 
A and B such that 



K 



(a,c) ■ 



(3) 

The initial parameters (a, c) can also play the role of being the final ones. In that case we use 
the parameters (a, c) , 

(4) 



It is a matter of substitution to verify that operators A and B intertwin the elements of K>(a,c) 
with those of JC^n and vice versa. Table 1 displays the set of basic intertwiners of diverse 
orders, all of them written in the confluent hyper geometric notation. 

Table 1: Diverse order intertwiners for the c.h.o. R x is the reflection operator R x <p(x) — tp(—x). 



Intcrtwincr 



Expression 



2 - c 



First 
First 
First 
First 
First 
First 
First 
First 

non-finite 
non-finite 



+c- 1 



(A- 1 ) 



a - 1 

a + 1 

a — c 

a — c + 1 

a — c -f- 1 

a - c + 2 



c+ 1 

c - 1 

c - 1 

c+ 1 

1 - c 

1 - c 

3 - c 

3 - c 



It is remarkable the presence of the zero order differential operator Q, which gives rise to 
a nontrivial intertwining operation: i( a -c+i,2-c) ( 5(a,c) = Q(a,c)Bu,c)- On the other hand, the 
non-finite order intertwiners V and W underlies some of the very basic relations obeyed by the 
c.h.f. such as the well known Rummer's first formula 0j . The first order intertwiners {A 1 , B l }f =1 
deserve as much attention because they factorize the c.h.o. £( a , c ) and £(a,c) : 

L( a ,c) = S (o,c)^(a,c) ~~ ^c)' ^(a,c) = ^{a,c)B\ a ,c) ~ Q(a,c)i i = 2 > 3 > 4 - ( 5 ) 

£(5,c) = ^( a ,c) B (S,g ~ ?(a,c)5 L (a,c) = B (a,c)^(a,c) ~ <?(a,c)> * = 1, 2, 3, 4. (6) 

where, by convention, we have implemented the induced action of the operators A and £> on 
R 2 as follows: A l (a,c) := A l ,Ja,c) = (5, c), B 1 (a,c) := £^ ac ^(a,c) := (a, c). The factorization 
constants <^ a c ^ are given by 

l\a,c)= a -^ <lfa,c)= a - 1 i <lia,c)= a - C > 9(o,c) = ° ~ L 

Now, in order to determine the explicit action of the free index operators on an arbitrary kernel 
fc( a ,c) > if is enough to consider a basis for K*( a ,c) as follows 

iFi(a,c;x), u(a, c; x) = x 1 ~ c \Fi (a — c + 1, 2 — c; x), (7) 
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where it is assumed that c ^ Z (the general case c £ Z is easily recovered afther an straight- 
forward limiting procedure [§]]). The action of the intertwiners Q, V, W and {A\ B l }f =l on this 
basis is displayed on Table 2. 



Table 2: Action of the basic free index intertwiners on the basis of K,i a 



Operator 


1 F\ (a, c; x) 


u(a, c; x) 


Q 


u(a + 1 — c, 2 — c; x) 


lFi(a + 1 - c, 2 - c; x) 



A 1 
B 1 
A 2 
B 2 
A 3 
B 3 
A 4 
B 4 



^) 1 F 1 {a, c+l;x) 
c — 1) \F\ (a, c — 1; x) 
c - 1) iFifa - 1, c- 
f) + 1, c + 1, x) 

a - c ) ii(a — c, 1 — c; 3:) 
^) u(a — c + 1, 1 — c; x) 
(c — l)u(a — c+1,3 — c, x) 
(c — 1) ii(a — c + 2, 3 — c; x) 



(1 — c) n(a, e+ 1; jc) 

( "i^-e 1 ) u ( a . c - li 

u ( a - 1, c - 1; 
(1 - c) u(a + 1, c+ 
(1 — c) i.Fi(a — c, 1 — c; 
(1 - c) i-Fi(a - c + 1, 1 - c;x) 
(f^) l F l< a - c+l,3-c;a) 
( "a-t 1 ) lJ ? l(i-c+2,3-c; a: ) 



V 



1 -Fi (c — a, c; 

— a, 2 — c;x) 



( — l) 1 c u(c — 0, c; x) 



As regards the composition of intertwiners, we first notice that Q 2 = 1. Hence, the kernel 
/C( a c ) is invariant under the twice iterated action of Q. The other operators are interrelated by 
means of Q: 

W = QV, A 2 = QA 1 Q, A 3 = QA 1 , A 4 = QA 2 , B 2 = QB 1 Q, B 3 = QB 2 , B 4 = QB 1 . (8) 

These expressions can be now used to compute the commutation rules obeyed by the intertwiners. 
On the other hand, according with Table 2, the functions i-Fi(a, c; x) and u(a, c; x) are such that 
they interchange roles under the action of Q. The same holds for other operators such that 
W or A 4 but adding a multiplicative non trivial constantf]. Thereby, Q becomes an important 
intertwiner in our approach, it works simply by transforming one element of the basis into the 
other one and vice versa. Such a behaviour is translated from the c.h.f. space into the space of 
parameters as a reflection operation. That statement is clear by considering a more convenient 
parametrization of R 2 , defined by the transformation 

a' = 2a-c, c ' = c-l, Va,cGR 2 . (9) 

With this new parametrization the operators Q and V are rewritten in such a way that their 
induced action on R 2 becomes linear, homogeneous and diagonal: 

Q(a',c') = (a',-c'), V(a', d) = (-a', J). (10) 

The induced action of the first order intertwiners in terms of the new labels is also simpler 

f A\a',d) = (a'-l,d + l), f A 2 (a',d) = (a'-l,d-l), 
\ -B 1 (a / , d) = (a' + 1, d — 1), \ B 2 (a', d) = (a' + 1, d + 1). 

This unpleasant situation can be solved by introducing maximal shape-invariant functions preserving their 
form under the action of all the intertwiners on Table 1 141. 
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A^a',d) = (a'-l,-d-l), f A\a',d) = (a' - 1, -d + 1), 
B 3 (a',c') = (a / + l,-c'-l), \ S 4 (a',c') = (a' + l,-c' + l). 

From 0) and Table 1 we see that V and VF are also reflection operators in the parameter 
domain. A different situation arises for the other intertwiners. For instance, according to 
the action of A 1 and B 1 on the point (a',c r ) € R 2 produces the displacements of a' and d . 
Therefore, by iterating the action of A 1 on R 2 we obtain (A ) n (a', d) = (a' — n,d + n), while 
for 5 1 it reads {B 1 ) m (a',c') = (a' + m,d — to). All these points are indeed in a straight line 
and form a linear discrete lattice on the space of parameters: {(a' + s,d — s),s G Z}. For each 
one of these points we can associate a kernel /C( a ' +s c /_ s ) and, in this way, the action of A 1 or 
B 1 on a c.h.f. /(a',c';x) can be understood as the mapping from /C( a / jC /) into /C( a /_ ljC / +1 ) or 
/C(a'+i,c'— l) respectively. The iteration procedure is now clear and the complete set of related 
c.h.f. {f(a r + s,d — s;x)}fL_ OQ is nothing but an invariant subspace under the action of A 1 
and Br. A similar situation occurs for the other first order intertwiners. In general, when the 
kernel of an operator X k £ {A k , B k }^ =1 is in a kernel /C( a ' iC '), the involved invariant subspace 
is bounded. We shall call an annihilation line of the operator X k , denoted by aZ[X fc ], to the set 
of points (a',c') G R 2 such that f(a',c';x) £ IC^^) iff X k a , c ,sf(a',d;x) = 0. From || and || 
it is easy to see that the factorization lines for the operators A k and B k are given by q k a , « = 
and q k a , = respectively. After a simple calculation and using (||) one gets the following 
relationships 

al[A 3 } = al[A\ al[A 4 ] = al[A 2 ], al[B 3 } = al[B 2 ], al[B A ] = al[B l ], 

which give the chance to construct a common set of functions T = {f(a' n ,c' m ;x)} representing 
an invariant space under the action of {A k , B k }f, =l . In this case the points {(a' n , c' m )} constitute 
a two dimensional discrete lattice on R 2 . The most interesting situation arises when there 
exists a critical point on the lattice because, by definition, this point is shared by two different 
annihilation lines. A situation giving rise to an invariant sector of bounded by that lines. In 
fact, there are four critical points on our (a',c') planef]: (—1,0), (0,1), (0,-1), and (1,0). We 
can then define the following sectors: 

Left Invariant Sector (L.I.S.). The point (a',c') = (—1,0) is at the intersection of aif-B 1 ] = 
{(a',c') | d = a' + 1} and al[B 2 ] = {(a',d) \ d = —a' — 1}. A left sector of the lattice is then 
bounded by these annihilation lines and it can be filled by successively applying higher powers 
of A 1 and A 2 on (-1,0). There is a left- corner c. hi. /(— 1,0; x) = 1 connected with this point 
while the related functions for the other points on the sector are of the form 

/(-l - n, n; x) = (-l) n , d > 0; /(-l - m, -m; x) = (-x) m , c < 0. 

In this way, we can characterize a support space generated by the basis functions {/(— 1 — m — 
n, m—n; x), m,n £ Z + }, called left invariant sector and denoted Tii, for which there is a (highest 
weight) irreducible representation of the algebra {A ,B = — (A 1 ) + , A 2 , B 2 = — (A 2 ) + ,Q} (for 
details see 0]). 

Right Invariant Sector (R.I.S.). The critical point (a',c') = (1,0) is at the intersection of 
allA 1 } = {(a', d) \ d = a' - 1} and al[A 2 } = {(a', d) \ d = -a' + 1}. The right sector bounded by 
2 The a' values run on the horizontal axis in a right hand frame of reference. 
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these lines is generated by means of B l ,B 2 . The involved c.h.f. are given by 



g(l + n, —n; x) = x n e x , c' < 0; 



g(l + m, m; x) = e x , c' > 0, 



from which, the right corner function g(l, 0; x) = e x can be obtained by making n = or m = 0. 
Therefore, we can construct a (lowest weight) irreducible representation of the algebra {A 1 , B 1 = 
— (A r ) + ,A 2 , B 2 = — (A 2 ) + ,Q}, whose support space spanned by {g(l+n+m, m— n; x), n, m G 
Z + } will be denoted Hr. 

Upper and Lower Sectors The critical point (0,-1) is at the intersection of aZ[A 3 ] and 
al[B 3 ] while (0,1) is at the crossing of aZL4 4 ] and al[B ]. However, there are not corner c.h.f. 
associated with that points and one is unable to construct another invariant sector (outside of 
the above mentioned ones) by means of any of the first order intertwiners A k or B k . Hence, 
there are only the two doubly bounded sectors already mentioned. 

It is now clear that only the left and the right sectors are of our interest. It is remarkable 
that the intertwiner Q becomes a self adjoint operator in the support spaces related with each 
of these sectors. Moreover, by taking a function (— l) n in the upper part of the left sector we 
get, after the action of Q, a function (— x) n in the lower part of that sector and vice versa. The 
same can be said about the functions living in the right sector. Hence, this operator intertwins 
the elements of a given sector with elements of the same sector. A very different situation holds 
for the reflection operator V which, by acting on a function (— x) m , gives e x x m . In other words, 
it maps the functions living in the upper part of the left sector into the functions of the upper 
part of the right sector and vice versa. The same situation occurs for the involved lower parts. 
The left and right invariant sectors are then intertwined by means of the operators Q and V 
whereas they are generated by the iterated action of the first order intertwiners on the left and 
right corner functions. 

At a first sight, it seems to be enough to consider the points on the lattice living in the left 
and right sectors. However, there is still a subset of points deserving attention. The straight line 
defined by the constraint d = —1/2 results invariant under the action of {A 3 ,B 3 }. This line cuts 
al[B 3 } on (-1/2,-1/2) and al[A 3 } on (1/2,-1/2). Hence, the points {(-1/2-n, -l/2)}~ =0 and 
{(1/2 + n, — 1/2)}^L will form the lattices for two irreducible representations of {A 3 ,B 3 }. A 
similar argument can be established for the line defined by c' = 1/2 and the operators {A 4 , B 4 }. 
Therefore, we have constrained the original (a 1 , c')-plane from all the R 2 into its more relevant 
part, composed by two invariant interrelated sectors plus two invariant lines. 

3 Schrodinger equations 

Let us start this section by introducing a free index operator M, written in the confluent hy- 
pergeometric notation as follows 



M f(a,c;x) := M M f(a,c;x), M^y-Y^o^ 1 (a,c;x) 



(13) 



where 



e: 



,x 



nl/2 



tp(a, c; x) = 



_x c y'(x)_ 
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and Y( a ,c) is an auxiliary operator changing the independent variable x into y 

x — >x(y), 

F(x) — ► = F(x(y)) . 



Yi 



(o,c) 



The action of the operator M on the kernel elements of £( a ,c) produces the mapping from 
the confluent hypergeometric equation Lu c \f(a,c;x) = to the time independent Schrodinger 
equation 



Mv) = Etj) 8 (y), 



(14) 



where 2m/T% 2 = 1 and the label s stands for the dependence on certain Schrodinger parameters. 
In other words, 

f(a,c)x) — ^—^My) = ^~ 1 ( a 5 c;x(y))/(a,c;x(y)). (15) 

It is easy to check that the identification of E and V(y) in (|i~4] ) depends exclusively on the 
specific analytical realization of y'(x), i.e., on the function we have taken as the new independent 
variable y(x). Once we know how to connect the confluent hypergeometric equations to a class 
of Schrodinger equations, we can translate all the studied mathematical properties into physical 
ones within the framework of the stationary Schrodinger wave- functions. Therefore, by using the 
first order intertwiners of Table 1 and M, we can establish the following commutative diagram 

M 



f(a,c;x) 



ips{y) 



>,c) 



B 



(S,c) 



(16) 



f(a,c;x) 



M 



My) 



where the operators A* and B| relate Schrodinger wave-functions with different labels: 

A* Mv) = MAJ 0)C) M~ l Mv) oc ^(y), 

b| My) = MBf a) gj m- 1 My) oc My)- 

Now, we proceed to identify the immediate physical potentials V{y) related with the c.h.f. 
by means of the above mentioned transformations. We shall focus on the main information 
drawing the readers attention to our previous work Q for details. In each case t stands for the 
angular quantum number I G Z + and E for the involved energy eigenvalues. Moreover, for all the 
following examples, the c.h.f. f(a', c'; x) are understood as appropriate linear combinations of the 
basis elements (|^) . Notice that it is always possible to choose between two different combinations 
of the pair (a', c'), which are labeled by an ± subindex in each case. In addition, we have clearly 
established the part of the (a' n , c^J-lattice where the involved c.h.f. live. The specific physically 
allowed values of these points can be obtained by asking for the square integrability property of 
the corresponding Schrodinger functions ip s . 



N-dimensional harmonic oscillator 



3) 



N > 2, 



a/2 
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I = d + + 1 - N/2 > 
E = -2a' + > 
upper part of the L.I.S. 

One-dimensional harmonic oscillator 



l=\-d_- N/2 > 
E = -2a'_ > 
lower part of the L.I.S. 



Vosdy) =y 2 , y = x 



1/2 



d + = 1/2 
E = -2a' + > 
upper left invariant line 



N-dimensional Coulomb potential 



{ c >_ = -1/2 
E = -2a'_ > 
lower left invariant line 



^cfoulM 



+ 



{21 + N- l)(2£ + JV-3) 

4y 2 



y cc x 



c'+2-N 



E = - 4- 



upper part of the L.I.S. 

Morse potential 



2-c' -AT 



lower part of the L.I.S. 



V&(!/) = (f) (e 2ay -2Ae^), y = (lnx) 1 /°, a > 0, A > 



f c' + = ^v^B > 1 



A = -a'_ 



+ 



upper part of the L.I.S. 



d =-l^E<i 

— a v 

A = -o'_ 

lower part of the L.I.S. 



As we can see, the same Schrodinger function ip s can be constructed by means of the c.h.f. 
parametrized by (a' + ,d + ) as well as that with parameters (a f _,d_). As d_ is always such that 
d_ = —d + , these functions live in different parts of the same (left) sector or line. Hence, it could 
be enough to consider either the upper or the lower part of that sector or line to describe the 
physical solutions of the related quantum problem. In other words, there is a mapping 2:1 from 
the kernel fC( a ,c) to the Hilbert space TC of square integrable functions ip s . A situation leading to 
a representation of Q in TL as an operator which does not transform the general solutions of the 
Schrodinger equation but changing only its parametrization. The same is true for the reflection 
operator V developing the first Kummer transformation on the c.h.f. 
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In summary, the action of the operator M on the c.h.f. (living only in the afore mentioned 
sectors) gives the stationary Schrodinger solutions for the oscillator, Coulomb, and Morse po- 
tentials. As we have seen, M connects the c.h.o. ([!]) with the Schrodinger equation related with 
these potentials. Hence, only the left and right invariant sectors, besides the invariant lines we 
are dealing with, are the responsible for bound states in the physical problems here considered 
(we will not discuss here the states corresponding to the continuum spectrum of these exam- 
ples.) This is why they have been called the physical sectors of the confluent hypergeometric 
functions space. On the other hand, the same mapping M allows the connection among the 
wave-functions of the diverese quantum systems we have presented. As a very last example, let 
us consider the N = 2 case for the oscillator and the Coulomb potentials. In order to fix the 
notation, we first write the energy of the discrete spectrum as follows: 

— 1 a 2 9 . . 

E = 2n , E C = (nc + 1/2)2 > E M = --u\ 1 (17) 

where O, C and M mean oscillator, Coulomb, and Morse, respectively. Thus, we have 
Oscillator-Morse: 

V^ M ' AM) (y) oce-^jT^W 2 ), n =\M,lo=VM- (18) 

Oscillator-Coulomb: 

^ c/c \y) oc y 1 /Vg l °'' o) (2y 1/2 /(2no+l) 1/2 ), n =2n c +l, t =2t c . (19) 

For other dimensions N of the oscillator or Coulomb potentials, the relations keep on valid, but 
with other restrictions on the values of (n,£). The above formulas (|i~8|)-([r9|) do not change the 
number of nodes but they change the local density probability so that there is a mathematical 
mapping among these three potentials which does not represent a physical equivalence. 

Acknowledgements: JN and LMN aknowledge financial support by DGES projects (PB94- 
1115 and PB98-0370) from Ministerio de Educacion y Cultura (Spain) and also by Junta de 
Castilla y Leon (C02/197 and C02/199). ORO acknowledges support by CONACyT project 
32086E (Mexico). 

References 

[1] S. Fliigge, Practical Quantum Mechanics (springer, Berlin, 1971) 
[2] G.A. Natanzon, Vestn. Leningr. Univ., Ser. 4: Fiz., Khim. 10 (1971) 

[3] A.F. Nikiforov and V.B. Uvarov, Special Functions of Mathematical Physics (Birkhauser, 
Basel, 1988) 

[4] J. Negro, L.M. Nieto and O. Rosas-Ortiz, J. Math. Phys. 41, 7964 (2000) 

[5] J. Negro, L.M. Nieto and O. Rosas-Ortiz, J. Phys. A: Math. Gen. 33, 7207 (2000) 



9 



